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ABSTRACT 

We present results of an analysis of the local {z ~ 0) morphology-environment relation for 911 
bright {Mb < —19) galaxies, based on matching classical RC3 morphologies with the SDSS-based 
group catalog of Yang et al., which includes halo mass estimates. This allows us to study how the 
relative fractions of spirals, lenticulars, and ellipticals depend on halo mass over a range of 10 11 ' 7 - 
10 14 ' 8 ft -1 Mq, from isolated single-galaxy halos to massive groups and low-mass clusters. We pay 
particular attention to how morphology relates to central vs. satellite status (where "central" galaxies 
are the most massive within their halo). The fraction of galaxies which are elliptical is a strong function 
of stellar mass; it is also a strong function of halo mass, but only for central galaxies. We interpret this 
as evidence for a scenario where elliptical galaxies are always formed, probably via mergers, as central 
galaxies within their halos, with satellite ellipticals being previously central galaxies accreted onto 
a larger halo. The overall fraction of galaxies which are SO increases strongly with halo mass, from 
~ 10% to ~ 70%. Here, too, we find striking differences between the central and satellite populations. 
20 ± 2% of central galaxies with stellar masses M + > 10 10 5 M Q are SO regardless of halo mass, but 
satellite SO galaxies are only found in massive (> 10 13 h~ x M ) halos, where they are 69 ± 4% of the 
M* > 10 10 ' 5 M satellite population. This suggests two channels for forming SO galaxies: one which 
operates for central galaxies, and another which transforms lower mass (M* < 10 11 M Q ) accreted 
spirals into satellite SO galaxies in massive halos. Analysis of finer morphological structure (bars and 
rings in disk galaxies) shows some trends with stellar mass, but none with halo mass; this is consistent 
with other recent studies which indicate that bars are not strongly influenced by galaxy environment. 
Radio sources in high-mass central galaxies are common, similarly so for elliptical and SO galaxies, 
with a frequency that increases with the halo mass. Emission-line AGN (mostly LINERs) are more 
common in SOs, but show no strong trends with environment. 

Subject headings: galaxies: active — galaxies: elliptical and lenticular, cD — galaxies: spiral — galax- 
ies: clusters: general — galaxies: groups: general — galaxies: evolution 



1. INTRODUCTION 



Galaxies come in many different shapes and sizes, 
but primarily consist of two dynamically stable compo- 
nents, bulges and disks, with additional contrib utions 
from quas i-stable features such as bars and rings. |Hub- 



ble ( 1926 ) devised what has become known as the "tun- 



mg fork" classification to describe this morphological 
schema: elliptical galaxies (~ pure bulge); spiral galax- 
ies (disks containing spiral features, both barred and un- 
barred, with a sequence Sa-S(j^]of decreasing bulge com- 
ponent and increasing spiral arm opening angle); and 
lenticular "SO" galaxies (defined by the presence of a disk 
with no discernable spiral arms). 

The abundance of these galaxy types is now known 
to correlate strongly with environment: elliptical galax- 
ies li ve preferential ly in regions of very high local den- 
sity (Dressier 1980), inhabiting the cores of clusters and 
groups rattier than their outskirt s (Mclnick fc Sargent 
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However, the total lummosi 



1991 
y-limi 



Wilman et al ||2009| 

ed traction ot ellipu 



cals is similar in a wide range of environments (Desai 
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evolved very weakly since z ~ 


1 (Dressier et al. 1997[ 


Fasano et al. 2000 Smith et al. 2005| Postman et al.| 



Lenticular galaxies, on the other hand, are globally 
no more abundant in clusters than they are in groups 
(Wilman et al. 2009), with a m uch weaker dependence 
on local density than ellipticals ( Dressier 1980 D ressier 
et al.||1997| |Postman et al.||2005|TFoggianti et al.||2008 ) 



A lower fraction at fixed luminosity is found only in the 
lower density field ( Wilman et al.|200"9 ). The global frac- 
tion of lenticular galaxies has grow n by a factor > 2.5 in 
groups and clusters since z ~ 0.5 dDr essler et al.||1997 



Fasano et al.|2000} |Wilman et al.|200P ), potentially more 



rapidly i n groupsand low mass cluste rs than in hig h-mass 
clusters ( |Poggianti et al. 2009| Just et al. 2010[ ). Since 
lenticulars and ellipticals make up the bulk of th e low 
redshift passive population (e.g. Bundy et al.|2010 ), this 
means that the majority of galaxies which have ceased 
forming stars since z ~ 0.5 have retained their disks as 
lenticulars. Beyond z ~ 0.5 (up to z ~ 1), there is no ev- 
idence for further evolution in the elliptical or SO fraction 
in the high-mass clusters s ampled to date (Smith et al. 
20051 i Postman et al.|2o05 1 . 



To build a global picture of galaxy evolution, it is im- 
portant to understand how the evolution of galaxies and 
their environment is intertwined, and which physical pro- 
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cesses drive the suppression of star formation and the 
morphological transformations which are required to ex- 
plain these observations. 

It is attractive to explain morphological evolution in 
the context of bulge growth through gal axy mergers 
(Springel et al. 2005 Hopkins et al. 2010), which take 
place in the center of halos as a natural consequenc e of 



hierarchical clustering and dynamical friction (e.g. De 
Lucia fc Blaizot||2007| |Pe Lucia et al.||2011| ). ElliptT 
cal galaxies almost certainly form through a mixture of 
major and minor mergers. However, the role of merg- 
ers in the formation of bulges in galaxies with disks is 
less clear. Cooling flows should form at the center of 
halos, and will reform a disk around the merger rem- 
nant so long as the gas is not too efficiently reheated. 
SPH simulations suggest that during mergers of galax- 
ies with high gas fractions, much of the gas retains its 
angular momentum and reforms a disk, and that some 
compone nt of the stellar disk can survive in all but 1:1 



mergers (Hopkins et al. 2009). Meanwhile, secular pro- 
cesses might co ntribute 5, 10% of the galaxy mass in a 
"pseudobulge" ( Kormendy fc Kennicutt|20 04), or poten- 
tially a larger traction at high redshitt wher e disks are 



less dynamically stable (Genzel et al. 2008). That SO 



spiral galaxies (e. 


2;. |Dressler 


1980 


Wilman et al. 


2009 


Laurikainen et al. 


2010) is consistent with a merger ori- 


gin. INonetheless, 1 


;here is currently no more than circum- 



stantial evidence that bulge growth and the suppression 
of star formation in most SOs are causally linked. 

Our observational picture still misses one vital and sur- 
prisingly straightforward ingredient: a full picture of how 
galaxy morphology depends upon environment at low 
redshift. Early studies of environment focussed almost 
exclusively on galaxy clusters, for which good statistics 
and strong trends are to be found even without high spec- 
troscopic completeness. The more recent emphasis has 
been on extending these trends to high redshift with high 
resolution Hubble Space Telescope (HST) imaging. Mod- 
ern local surveys such as the Sloan Digital Sky Survey 
(SDSS) are so large that much effort has been invested 
in the automated classification of morphological prop- 
erties, such as the simple SDSS "concentration" index, 
and multi-component decomposition with diff erent levels 



of sophistication (re solution dependent, e.g. Allen et al. 
2006 Gadotti|2009[ ). An interesting alternative solution 
is the "Galaxy Zoo , providing visual classifications for 
the whole SDSS, by employing an en thusiastic public to 
classify galaxies (Lintott et al. 20111. 

The disadvantage of all these approaches is that new 
classification schemes tend to be relatively simplistic, 
with unknown systematics, and produce data which is 
difficult to assess on the basis of our existing understand- 
ing of the Hubble scheme. Standard Hubble type classi- 
ficat ions within SPSS are available only for subsamplcs 
(e.g. [Fukugita ct al. 200 71), alth ough these are becoming 
larger qiNair fc Abraham||2010a| ) . 

In this paper we have taken the simple step of tak- 
ing a sample of galaxies fro m the Third Reference Cata- 



log of Bright Galaxies (RC3, de Vaucouleurs et al.|1991 1 
which is matched to the SPSS Data Release 4 (DR4) to 
provide classical morphological information for a large 
sample of galaxies for which the selection function and 
environment can be properly characterized. For our pur- 



poses this provides an ideal sample to study environmen- 
tal trends for galaxies with detailed and well understood 
classifications. 

In section [2] we describe the sample, discussing the se- 
lection function of RC3 galaxies and the group catalog 
which we use to describe galaxy environment. Section [3] 
presents our results, including Hubble type fraction as a 
function of stellar and halo mass, a dichotomy between 
central and satellite galaxies, the nature of activity in the 
SDSS fiber spectra of galaxies with different morphologi- 
cal types and environments, and the dependence (or lack 
thereof) of morphological components (bars, rings etc) 
on environment. This is further discussed in the context 
of a hierarchically evolving Universe in section [4j and in 
section [5] we present our conclusions. 

To compute distances, absolute magnitudes, and stel- 
lar masses, we assume a ACDM cosmology with Q m, A, 
and HO equal to 0.3, 0.7, and 71 km s _1 Mpc -1 , respec- 
tively. Halo mass es are presented in /i _1 M , to retain 
the native units of Yang et al. ( 2007| )P| 



2. SAMPLE 
2.1. The RC3 Catalog 

The RC3 catalog (Third Reference Catalog of Bright 
Galaxies jj provi des information for a large sample of 
nearby galaxies (de Vaucouleurs et al. 1991). It "at- 
tempts to be reasonably complete for galaxies having ap- 
parent diameters larger than 1 arcmin at the D25 isopho- 
tal level and total 5-band magnitudes Bt brighter than 
about 15.5, with a redshift not in excess of 15,000 km/s". 
Some smaller, fainter, and more distant galaxies are also 
included. The most important aspect from our point of 
view is that detailed morphological classifications (Hub- 
ble types) from photographic plates are provided. 

Inspection of SDSS imaging for a subsample of RC3 
galaxies confirms the majority of classifications. How- 
ever, a significant minority of galaxies appear to be 
falsely classified, usually in the sense that some galax- 
ies appearing to be early-type spirals or ellipticals are 
instead classified SO. We describe our reclassification 
process in section |2.4| 



2.2. The SDSS Catalog 

The S PSS DR4 (Sloan Digital Sky S urvey Data Re- 
lease 4, Adelman-McCarthy et al.|2006 1 provides ugriz 
photometry and spectroscopy for 565,715 galaxies across 
a total of 4783 square degrees. In the main sample, this 
is highly complete down to a limiting magnitude of 17.77 
in r-band dereddened Petrosian magnitudes, and limit- 
ing surface brightness of /i r = 23.0 mag arcsec -2 . Low 
level incompleteness exists mainly in the highest density 
regions due to the inability to assign fibers to targets 
with separations < 55". 

The SPSS spectroscopic sample is incomplete at the 
bright end due to fiber magnitude limits, applied to 
avoid saturation and excessive cross-talk in the spec- 
trographs. Our catalog is based on the PR4 version of 
the New Yor k University Value A dded Galaxy Catalog 
(NYU-VAGC iBlanton et"aT1[2005l )), which includes ob- 
jects with bright fiber magnitudes and is photometrically 

4 h = 0.71 in our adopted cosmology. 

5 http:/ /heasarc. nasa.gov/W3Browse/all/rc3. html 
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recalibrated across the sky. If we restrict ourselves to the 
SDSS area with DR4 spectroscopic coverage, we measure 
redshift completenesses of ~ 93% for r < 16, ~ 88% for 
13 < r < 14, and - 80% for 12 < r < 13 (i.e., the frac- 
tion of galaxies in the NYU-VAGC which a ctually have 
SDSS spectroscopy). Strauss et al. (2002) studied the 
incompleteness of SDSS galaxy spectroscopy using local 
catalogs and found typical incompleteness of ~ 5% for 
bright (r < 15) galaxies. Significantly, this incomplete- 
ness was due to overlaps with saturated stars. What this 
indicates is that the bright-galaxy incompleteness is both 
low and, crucially, not related to galaxy type. Therefore 
this incompleteness is not important for our analysis. 

2.3. A Matched RC3-SDSS Catalog 

We use the NYU- VAGC match to RC3 galaxies from 
Blanton et al. ( 2005 ) . Although the matching radius was 
45" to compensate for the variable RC3 astrome trical ac- 
curacy, all matches are within 15" (figure 7 of Blanton 
et al. 2005). Most of the galaxies under consideration are 
themselves larger than 15" in size, and thus the incidence 
of false matches is low. There may be few occasions in 
which a false match is made for interacting or neighbour- 
ing galaxies, but this is unlikely to greatly influence our 
overall statistics. 

Our full matched RC3-SDSS catalog contains 1340 
galaxies, mostly in the range 0.01 < z < 0.04 (median 
z = 0.023), with a tail of bright galaxies extending to 
2 = 0.13. 

2.4. Morphological Reclassification 

Subsequent examination of individual galaxies sug- 
gested that some of the RC3 classifications were in error 
(e.g., a "spiral" galaxy that was clearly an elliptical, or 
an "SO" that had strong spiral arms). To address this, 
we carried out a reclassification analysis for the sample. 
Two student interns independently examined all galaxies 
in our sample with assigned halo masses, identifying 406 
cases with potentially incorrect general classification (i.e. 
E, SO or spiral). We then examined each of these objects 
in detail, reclassifying them if necessary)^] In cases where 
a lenticular was reclassified as a spiral or vice-versa, we 
carried over any original disk-structure notations (bars 
and rings) from the original RC3 classification. 165 of 
the 406 flagged galaxies ended up with new classifica- 
tions, including 6 of 23 ellipticals, 80 of 129 SOs and 65 
of 230 spirals (a few galaxies have irregular or merger 
types). One galaxy (UGC5677) was clearly matched to 
the wrong SDSS object, due to especially poor RC3 as- 
trometry; we assigned it zero weight in order to remove 
it from the sample. Finally, we also classified a total of 
55 galaxies wit h B < 16 (our primary magnitude limit; 
see Section 2.5) which had no classifications in the RC3. 

The most striking aspect of this analysis was the high 
fraction of reclassified lenticulars: a total of 40% of 
the original (RC3) lenticulars ended up with different 
Hubble types (16 as ellipticals, 34 as spirals, and 9 as 
peculiar/merger-remnant systems). A similar high pro- 
portion of reclassified RC3 lent iculars can be see n in Fig- 
ure 1 of |Fukugita et al.| (|2007[) and Figure 14 of |Nair fc 



Abraham (2010a). This led to a general reduction in the 
overall traction of SOs and smaller increases in the el- 
liptical and spiral fractions. The reclassification process 
actually strengthened the trends described in section [3] 

2.5. SDSS Photometry and the RC3 Selection Function 

Our goal is to examine the morphological composition 
of galaxies using the corrected RC3 classifications. This 
requires a detailed understanding of the inherently non- 
uniform RC3 selection, using the full SDSS catalog as 
the nominally complete "parent sample" . 

Our photometry is based on the model magnitudes in 
the NYU-VAGC c atalog, which we h ave corrected for 
galactic extinction (Schlegel et al. 19981. For large galax- 
ies (r 50 > 10", where r 50 is the radius containing half 
the r-band Petrosian flux) , the SDSS background can be 
oversubtracted, due to the outer part of the galaxy being 
treated as part of the sky; this effect is approximately 
color independent (Blanton et al. 2011). To account for 
this, we applied the photometric correction published by 
Blanton et al., as a function of the measured half-light 
radius r 50 (their Table 1, correction to v5.4 sky). This 
cor rection is g e nerall y co nsistent with those publi shed 



by | West et al.| ( |20i0j ) and |Hyde fc Bernardi| p009| . A 
comparis on with" the latter cor rection tor ellipticals - Fig- 
ure 14 of Blanton et al. (2011 ) - suggests that these cor- 
rections might depend upon morphology. However, the 
differences are only significant (d(mag) ~ 0.5) for galax- 
ies with rso > 25". Since only one of the galaxies in 
our sample is this large, we do not consider morphology- 
dependent corrections to be something to worry about. 

As noted above, the main selection functions apply- 
ing to the RC3 catalog are limits on (photographic) blue 
magnitude and isophotal diameter. To determine how 
the magnitude limit affected the selection, we generated 
an equivalent B magnitude for all SDSS galaxies, using 
the SDSS u and g magnitude^ 



B 



0.8116(u - g) + 0.1313. 



(1) 



° For reclassification, both color JPEG and background- 
subtractcd p-band FITS images were examined. 



To evaluate the effects of the diameter limit, we used the 
SDSS rgo measurement (the radius containing 90% of the 
Petrosian flux in r-band). 

How well do our SDSS-derived B magnitudes match 
the (largely photographic) magnitudes used as limits for 
RC3? Measurements of the RC3 total -B-band magni- 
tude, Bt, are available for a sub-sample of 150 SDSS- 
RC3 galaxies. For these objects we compare Bt to B 
derived from the SDSS photometry. We find that the 
RC3 and SDSS b magnitudes are comparable in 144/150 
cases, with significant scatter due to the large RC3 pho- 
tometric errors. The other 6 objects have significantly 
underestimated fluxes from SDSS, with measured g-band 
magnitudes of g > 15.8 (corresponding to B > 16.0). 
Five of these are clearly cases where deblending-related 
problems have affected the SDSS measurements. The 
sixth is an RC3 object (NGC842) mis-matched to the 
wrong SDSS galaxy due to a ~ 1' error in its RC3 cata- 
loged position. We remove these objects by only consid- 
ering galaxies with g < 15.8 and B < 16.0, although the 
"true" magnitudes might be brighter than this limit. As- 
suming that this analysis is representative, we estimate 

7 http:/ /www. sdss.org/dr7/algorithms/sdssUBVRITransform. html 
#Lupton2005 
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that ~ 4% of objects are lost due to poor SDSS photom- 
etry, and a < 1% occurance of RC3-SDSS mismatches. 

The first two panels of Figure [I] illustrates how SDSS- 
RC3 galaxies and the (almost complete) SDSS parent 
sample populate the magnitude-size parameter space. 
Not surprisingly, fainter or smaller galaxies are less likely 
to be found in the RC3 catalog. The third panel shows 
the ratio by number of RC3 to parent-sample objects as 
a function of these two parameters. To correct for the 
RC3 selection biases, we weight SDSS-RC3 galaxies by 
the inverse of this selection function: weight(B, rg ) = 
Npare n t(B ,r 9 o) com p U t ec i m \,[ ns f o.5 mag in B mag- 

nitude and 5" in rgn. Some of the fainter, smaller ob- 
jects have weights of up to 41. A threshold brighter than 
B = 16.0 would remove these objects, but significantly 
worsen the overall statistics. We examine the robust- 
ness of our results by applying three different limits - 
B < 16.0, 15.5, and 15.0 - and checking that each result 
is consistent with each of the three limits applied. 

We would like to compute fractions of the complete 
population down to a given luminosity limit. However, 
more luminous galaxies are visible to larger distances, 
and thus over larger cosmological volumes (this is the 
well-known Malmquist bias). To correct for this, we com- 
pute V/V max weights for our galaxies. In practice, V max 
is the volume within which a given galaxy is visible, and 
V is the volume of the survey, defined to be the maximum 
value of V max (minimum V/V max = 1). To ensure we are 
insensitive to large V/V ma , x weights, we have carefully ex- 
amined our results both with and without these weights 
applied, and as a function of luminosity. Robust results 
are achieved with a luminosity limit of Mb — —19.0. 

2.6. The Group Catalog 

For a full picture of how galaxy morphologies de- 
pend upon e nviro nment, we use the group catalog of 



Yang et all ( 20071 Y 07) constructed from SDSS DR4, 
and applied to galaxies in the range 0.01 < z < 0.2. 
The Y07 catalog is constructed iteratively using an al- 
gorithm which is an amalgam of friends-of-friends link- 
ing of galaxies to form groups and abundance matching 
to assign masses to those groups. Group halo masses 
are assigned based upon the rank order in terms of the 
group total stellar mass or luminosity of all galaxies 
brighter than an evolution and k-corrected r-band ab- 
solute magnitude of —19.5. We use the stellar-mass- 
based halo masses, for which halos are complete down 
to Mhaio = 10 11,63 h _1 M ; this corresponds to a single 
red galaxy of stellar mass M* = 10 10 M just making 
the —19.5 cut in luminosity. This method provides halo 
masses for single as well as grouped galaxies. 

We use Sample II as described by Y07, which contains 
a total of 369,447 galaxies including those with redshifts 
from sources other than SDSS and assumed redshifts for 
fiber collision galaxies as described. Galaxies and groups 
with incomplete local information are removed from the 
catalog, which means galaxies with comple teness in the 
local "tiling sector" ( iBlanton et al.||2005[ ) of C < 0.7 
(missing > 30% of neighbours) and groups on the edge 
of the survey with completeness f e d ge < 0.6 (> 40% of 
group members with C < 0.7). Isolated galaxies with 
stellar masses M* < 10 10 M are not assigned halo 



halo mass limit. The vast majority of such galaxies are 
fainter than our Mb = —19.0 cut. 

Figure [2] examines how halo mass estimates for groups 
with > 3 members correlate with line-of-sight velocity 
dispersion a\ os (top) and number of confirmed members 
Njnem (bottom), dies is computed using the "Gapper" 
algorithm, a ppropriate for gro ups with a small number 
of members ( Beers et al.|1990 ). The overplotted median 
Mhaio-cios relations (using a running bin of 50 galaxies - 
or 5 galaxies for the most massive 25 groups) shows that 
a typical group halo mass of ~ 10 13 hr 1 M corresponds 
to a measured velocity dispersion of a\ os ~ 135 kms -1 
and iV mem ~ 3-18 (median ~ 6). These masses are 
~ 40% less than a simple prediction derived from the 
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However, this is expected; ap- 
plication of the group-finder algorithm to mock cata- 
logs predicts this ~ 40% discrepancy which results from 
the trimming of high peculiar velocity members at the 
tails of the velocity histogram (Y07). Thus Mhaio ~ 
10 13 h _1 Mq corresponds to a true velocity dispersion of 
<7ios ~ 190 kms -1 , and galaxies with large peculiar ve- 
locities will typically be assigned their own halo. 

Y07 also classify galaxies as either central or satellite. 
In reality, a "central" galaxy is just the most massive 
galaxy in any given halo (i.e. equivalent to the Brightest 
Group Galaxy, BGG), regardless of its actual physical 
location within the gro up, and all other gal axies in the 
group are "satellites". |Skibba et al. (20111 show that 
a high fraction of so-called "central" galaxies are not in 
fact located at the center of the group in terms of ei- 
ther projected position or velocity: going from ~ 25% 
for IO^H^Mq < M h aio < lO 13 h- 1 M to - 40% for 
Mhaio <; lO 137 h _1 M . This must be a consequence of 
the dynamical state of many of these systems which will 
be caught in the midst of a major halo merger, or at a 
time before relaxation of the descendent halo. Nonethe- 
less, there are good theoretical and observational rea- 
sons to treat the most massive galaxy differently from 
the other galaxies: it should have the deepest potential 
well of any galaxy in the halo, especially if it does sit at 
the center of the global potential (when one exists after 
relaxation). This means halo gas can cool more easily 
onto this galaxy than it can onto any of the others. Ob- 
servationally, these two types of galaxy have been shown 
to behave quite differently on average: central galaxies 
are more likely than satellites to be forming stars and/or 
hosting radio AGN at fixed stellar and halo mass flWein- 
|mann et al.||2006| 



von der Linden et al. 1120071 [Best et al. 



masses, since they will live in halos below the 10 



11.635 



M M 



j]von der Linden et al.|2010[ ) . Therefore, we use these 
classifications to examine the dependence of galaxy mor- 
phology on whether the galaxy is a central or satellite 
galaxy as well as on the underlying halo mass. 

Of the 1340 galaxies in our SDSS-RC3 sample, 1194 
have B < 16; the sample size shrinks to 1064 when our 
Mb < — 19 cut is also applied. Within this subsample, 
there are 911 galaxies which have halo masses assigned by 
Y07, with 729 of these being central galaxies and 182 be- 
ing satellites. It is important to note that our SDSS-RC3 
sample does not include the richest, most massive galaxy 
clusters: the maximum halo mass is ~ 7 x 10 M , while 
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Fig. 1. — 2D histograms used to determine the RC3 selection function, using a B = 16 limit. The left and middle panels show counts 
of RC3 galaxies and SDSS parent sample galaxies binned by rgo (the radius containing 90% of the r-band Petrosian flux) and B-band 
magnitude; horizontal color bars indicate the scales in galaxies per bin. The right panel shows the ratio of the left and middle panels, 
which is the basis of the RC3 selection function. The horizontal color bar indicates the ratio (from 0.0 to 1.0); white bins indicate no data 
(i.e., no galaxies in either sample are found in those bins). 

galaxies tend to have lower B-band stellar M/L light ra- 
tios, although the highest M/L ratio galaxies are also spi- 
rals (with dust-obscured B-band light). Ellipticals and 
SOs populate a fairly tight relation, indicating that there 
is little difference between a luminosity-based cut and a 
mass-based cut for E versus SO classifications. At our 
selected luminosity cut of Mb = —19.0, an early-type 
galaxy has a typical stellar mass of M* ~ 10 1 M©. 
Cutting in mass instead of luminosity would reduce the 
unweighted total spiral fraction by less than 4%. 

Figure [3] also shows that spiral galaxies with Mb < 
— 19 can have masses down to ~ 10 9 5 M Q . Because the 
halo mass estimates are ultimately based on galaxy stel- 
lar masses, these different stellar-mass limits have impli- 
cations for our ability to find different galaxy types as a 
function of halo mass and galaxy status (central versus 
satellite). The simplest case is for central galaxies. If 
we consider isolated galaxies (where the galaxy is auto- 
matically the central galaxy of its halo), then a stellar- 
mass limit of 10 10 " 5 M© implies a halo mass limit of ~ 
10 11 ' 8 hr 1 M© for central E/S0 galaxies. The lower mass 
limit for spiral galaxies means that they can be found as 
central galaxies for halos with lower masses, down to the 
Y07 completeness limit (M halo = 10 11 - 63 h" 1 M ). This 
means that for halos masses < 10 12 M©, we should 
expect to find few or no central E/S0 galaxies, solely due 
to our luminosity cut. 

A similar, albeit more complicated, effect applies to 
satellite galaxies. The minimum halo mass for a satel- 
lite galaxy can be estimated assuming that it is the sec- 
ond most massive galaxy in a two-galaxy group, with 
the central galaxy only marginally more massive. In 
this case, a 10 1 M© satellite galaxy will have a halo 
mass > lO 12 5 hr 1 M©, and we should expect to under- 
count E/S0 satellite galaxies in less massive halos. A 
spiral satellite with M* ~ 10 95 M©, on the other hand, 
could reside in a halo with total stellar mass content 
M* > 10 9 ' 8 M© corresponding to halos of low mass (po- 
tentially down to the Y07 completeness limit). 

The fact that our Mb < — 19 cut largely excludes E 
and SO galaxies with stellar masses < 10 10 5 M©, while 
still including a large number of spiral galaxies with 
smaller stellar masses, means that direct comparisons of 
E, SO, and spiral galaxies will be biased for stellar masses 



the maximum velocity dispersion is ~ 740 km s 1 . 

2.7. Stellar Masses 

To examine how the galaxy population varies with stel- 
lar mass, we calculated stellar masses for each gal axy us- 
ing the color -based mass-to-light (M /L) ratios of Zibetti 
et al.| ( |2009[ hereafter Z09), using SDSS g - i colors and 
i- band absolute magnitudes (including the necessary k- 
corrections) . We prefer this app roach over using the stel- 
lar masses of |Yang et all ( |2007D which are based on the 
M/L ratios ot |rjell et la!] ( |2003p , because the Zibetti et 
al. M/L ratios include the effects of dust reddening and 
extinction, along with a spread of possible star-formation 
histories. 

We also investigated u sing the stellar mass derivations 
of Gallazzi et al.| ( |2"005| hereafter G05)(^] which take ad- 
vantage ot SUSS spectroscopy for each galaxy. The G05 
masses were estimated by fitting five spectral absorption 
features to estimate a z-band M/L ratio, applied to the 
z-band luminosity. The drawback of this approach, from 
our perspective, is that the SDSS fiber spectra sample 
a relatively small, central region of the galaxies: since 
RC3 galaxies usually have diameters > 1', the 3-arcsec 
SDSS fiber aperture captures only the inner region of the 
galaxy. For spiral galaxies, this can mean that the spec- 
troscopic M/L reflects the bulge or nuclear region, rather 
than the galaxy as a whole. In addition, spectroscopy- 
based stellar masses are not available for some galaxies 
(142 of the 1340 galaxies in our full RC3-SDSS matched 
sample) . 

In the rest of this paper, we mainly concentrate on the 
Z09-based stellar masses, but we tested all trends with 
stellar mass against the G05 and Y07 masses as well, and 
note when these provided different results. 

2.8. Systematics and Limits in Stellar and Halo Masses 

Figure [3] illustrates how galaxies of different morpho- 
logical type populate the M^-M* plane. (The B-band 
luminosity Mb is computed applying both distance mod- 
ulii in our chosen cosmology and k-corrections based on 
the kcorrect code of Blanton & Roweis 2007.) Spiral 

8 These stellar masses, as well as Ho linewidth measurements, 
were taken from the MPA public webpage http://www.mpa- 
garching.mpg.de/SDSS/DR4/ 
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Fig. 2. — Y07-calibrated group halo mass versus group rest-frame velocity dispersion (top) and number of group members (below) for 
groups with at least 3 members. The red line indicates the median of the distribution as a function of halo mass, computed over a running 
bin of 50 galaxies (5 for the 25 highest mass groups, dotted line). The blue solid line in the top plot is the result of a simple virial prediction, 
and the magenta dashed line is the simple fitting function used by Y07. 

spectroscopic properties, and morphological substructure 
(bars and rings) - depend on global galaxy properties and 
on group properties. 



< 10 M . Consequently, when we consider fractions 
of different morphological ty pes a s a func tion of stellar 
or halo mass (e.g., in Sections 3.2 and 3.4 1, we limit our- 
selves to M* > 1O 1O 5 M . This limit does not apply 
when we consider fraction within a given morphological 
class - e.g., the fraction of ellipticals with optical AGN 
spectra, or fractions of spirals with outer rings. 

3. RESULTS 

In this section we look at how various galaxy classifica- 
tions - broad Hubble types (elliptical, lenticular, spiral) , 



3.1. Uncertainties and Significance Testing 

We use fractions to discuss the morphological com- 
position of the local galaxy population. Fractions are 
by definition relative quantities and as such allow com- 
parisons between different types, doing away with the 
need to renormalize to the total galaxy population in 
a given bin. However, the reader should keep in mind 
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Fig. 3. — Stellar mass (M*) versus B-band luminosity (Mj) for all galaxies in our sample. Elliptical, 
by red circles, green triangles, and blue stars, respectively. The vertical dashed line marks our Mg = 



SO, and spiral galaxies are indicated 
— 19 luminosity cutoff. 



that focusing on fractions does have its own problems. 
For example, the total galaxy population in a given bin 
of stellar or halo mass, or luminosity, is not necessarily 
conserved with time. Since we consider multiple mor- 
phological types (e.g., elliptical versus SO versus spiral), 
the fraction of a given type can change due to transfor- 
mations between other types, and not just because the 
total number density of the first type is increasing or de- 
creasing (e.g., if the SO number density stays constant 
but spirals merge to form ellipticals, then the total num- 
ber of E+Sp galaxies decreases and the SO fraction will 
increase). 

Our plots show fractions for each specified type within 
various bins. The err or bars in these plots are 68% confi- 
dence limits from the Wilson (1927) binomial confidence 



galaxy is SO or not, galaxy is barred or not) as a func- 
tion of some independent variable x is modeled using the 
following function: 



P 



1 



1 + e - a+ P x 



(2) 



interval, which is a more accurate way of estimating bi- 
nomial uncertainties than the commonly used Gaussian 
approxim ation, especially as the frequencies approach 
or 1; see Brown et al. (2001 1 for a comprehensive discus- 
sion of binomial uncertainties. 

In the case of weighted counts, there is less guidance 
on the proper way to compute uncertainties. We esti- 
mate the uncertainties by first rescaling all (weighted) 
counts so that the total counts are equal to the original 
(unweighted) total counts in a given bin, and then com- 
puting the Wilson confidence limits using the rescaled 
counts. 

We estimate the significance of apparent trends in 
these plots via a logistic regression analysis. In logistic 
regression, the probability of a binomial property (e.g., 



where P is the probability for a galaxy having the par- 
ticular property. The coefficient a corresponds approx- 
imately to an intercept value, while f3 is analogous to a 
slope and measures how strong the trend is. Note that 
the function always has values between and 1 (appro- 
priate for a probability), and is either monotonically in- 
creasing or decreasing. 

In principle, we could fit straight lines to the binned 
frequencies and use the null probability for a nonzero 
slope as an estimate of a trend's significance, but the lo- 
gistic approach has two key advantages. First, a linear fit 
will yield (meaningless) frequency values above 1 or be- 
low at some point, while the logistic curve is bounded 
between and 1. Second and perhaps more importantly, 
linear fits to frequencies will be biased by the specific 
binning scheme used (number and spacing of bins) , and 
by the assumption of Gaussian errors - an assumption 
that is not necessarily true for frequencies, especially for 
frequencies close to or 1. Logistic regression uses all 
data points individually, making binning and error as- 
sumptions irrelevant. 

The specific logistic regression code we use is from the 
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Surv ey p ackager) implemented in the R statistical lan- 
guage^] this package allows for individual data points to 
have weights. The fitting process yields null-hypothesis 
probabilities for the intercept a and slope ft (i.e., the 
probability that a logistic model with a value of for the 
given coefficient could explain the data); we use the null- 
hypothesis probability for (3 (Pp=o) as an estimate of an 
apparent trend's significance. We also quote the best- 
fit value of ft and its uncertainty; positive values mean 
the frequency increases with the independent variable, 
while negative values indicate a trend of decreasing fre- 
quency. We caution the reader that this is not a univer- 
sal indicator for a trend: trends more complicated than 
monotonic, smooth increase or decrease in probability 
could be poorly fit by the logistic model. Nonetheless, 
we feel this is superior to the commonly used approach of 
least squares linear fitting of binned data, for the reasons 
given above. A two-sample test such as the Kolmogorov- 
Smirnov test is also not advisable, both because it is 
calibrated only for unweighted data and because it mea- 
sures the maximum, local difference between two samples 
(e.g. ellipticals and non-ellipticals) rather than the sig- 
nificance of a global trend. 

3.2. Morphological Fractions Versus Luminosity and 
Stellar Mass 

Figure [4] shows how the fractions of elliptical, lentic- 
ular, and spiral type galaxies depends upon _B-band lu- 
minosity within the full SDSS-RC3 sample. Selection 
weights are applied to galaxies; when plotting against 
luminosity, it is not necessary to apply y/Vmax weights. 

The overall elliptical fraction is low, but clearly in- 
creases with luminosity ((3 — — 1 ± 0.24, i-^o = 1.3 x 
10 -5 ; note that this indicates a frequency that decreases 
as Mb becomes more positive). The observed fraction 
of SO galaxies decreases with luminosity (/? = 0.8 ± 
0.3,Ps = o = 0.0088), whilst the fraction of spiral galax- 
ies show little, if any dependence on luminosity down to 

M B 19 (ft = -0.12 ± 0.22,P^ =0 = 0.59). From a 

theoretical standpoint, this is perhaps surprising, indi- 
cating that the overall probability that spiral arms have 
faded (likely associated with the suppression of star for- 
mation in a disk galaxy) is not a strong function of its 
luminosity. 

We examine the stellar mass dependence of morpho- 
logical fractions in the top panel of Figure [5j cut at 
Mb = —19 and weighted to account for both selec- 
tion and volume (Vmax)- We also restrict our plots, 
and our logistic regression analysis, to galaxies with 
M* > 1O 1O ^M , because selection effects can lead to 
artificial suppression of elliptical and SO popul atio ns rel- 
ative to spirals below this mass limit (Section |2.8[ ). The 
fraction of ellipticals does increase significantly with in- 
creasing stellar mass [p = 3.2 ±0.64, P p=0 = 4.3 x 10~ 7 ). 
Overall, the SO fraction appears to decrease as stellar 
mass increases, although this is not formally significant 
09 = -0.32 ± 0.74, P 0=o = 0.66). The fraction of spiral 
galaxies is essentially constant at ~ 65% over this mass 
range (/3 = -0.74 ± 0.68, Pp= = 0.27), similar to what 
was seen for trends with Mb- 

9 http:/ /faculty.washington.edu/tlumley/survey/ 

10 http://www.r-project.org/ 



3.3. Morphological Fractions Versus Halo Mass 

We present the fraction of galaxies with Mb < —19.0 
and > 10 10 5 M as a function of group/halo prop- 
erties in Figure [6j The top panel presents the trends 
as a function of halo mass. This shows a strong depen- 
dence of morphological type on halo mass. The elliptical 
fraction is practically zero below Mhaio ~ lO 12 h _1 M 
(but this is most like ly due to the selection effect dis- 
cussed in section 2.8), then rises to a roughly constant 
level of ~ 5-10% tor more massive halos. The lenticu- 
lar fraction increases dramatically from ~ 10% in the 
very lowest halo masses to ~ 80% at the high-mass 
end (the dashed green line shows the logistic fit, which 
has /? = 1.1 ± 0.32, P/3 =0 = 0.00071. The fraction of 
spiral galaxies naturally compensates for these trends 
with halo mass, decreasing from ~ 85% to ~ 15% go- 
ing from the lowest mass to the highest mass halos 
08 = -1.2 ± 0.31, P/^o = 0.00013 for the logistic fit 
shown by the dashed blue line). Significant trends for 
SOs and spirals are also found with group velocity dis- 
persion <7i os and number of members N mem as the inde- 
pendent variable (logistic fits, shown in the center and 
lower panels). 

The most dramatic outlying point in Figure [6 is the 
highest bin of N mcm (bottom panel) , in which the spi- 
ral fraction jumps to almost 80%. This bin is populated 
by a single group , the Abell 2199 supercluster (see e.g. 
Rines et al. 2002 1, wh ich has apparently been merged by 
the Yang et al. 1 2007 ) group-finder algorithm into a sin- 



gle, overmassive group with 433 members. Most of the 18 
RC3 galaxies belonging to this system lie outside the col- 
lapsed regions of the supercluster (including Abell 2197 
and Abell 2199). Their spiral morphologies are therefore 
not surprising. 

3.4. Comparing Central and Satellite Galaxies 

Figure [7] splits the galaxy population into central and 
satellite galaxies, which (at least theoretically) should 
be subject to different physical processes. The mor- 
phological fractions are plotted against both the stel- 
lar mass of individual galaxies ( M* , left panels) and 
halo mass ( Mhaio j right panels). Once again, we restrict 
the plots and logistic-regressi on a nalyses to galaxies with 
M* > 10 10 5 M m (see Section"" 



3.2| 



The dependence of galaxy morphology on galaxy stel- 
lar mass shows some interesting differences when we 
consider central and satellite galaxies separately. The 
roughly constant fraction of spirals as a function of stel- 
lar mass for all galaxies (Figure [5]) is replicated by the 
central galaxies; but for satellites we evidence of a trend 
where the frequency of spirals increases with mass, at 
least up to M* — 1O 115 M . For SOs, the apparent 
differences are even stronger: central SOs are a roughly 
constant fraction at all masses < 1O 115 M , while satel- 
lite SOs show a steep drop in frequency as stellar mass 
increases (however, this trend is not formally signfi- 
cant: j3 = —3 ± 1.8,P^ = o = 0.1). Clear trends are 
harder to discern for ellipticals; in general, they resem- 
ble the elliptical trend for all galaxies (Figure pi), with 
frequency increasing weakly but significantly with stel- 
lar mass (ft = 3.1 ± 0.64^=0 = 1.2 x 10~ 6 for central 
ellipticals) . 

When we turn to the question of how the morphological 
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Fig. 4. — Fractions of elliptical (red circles), lenticular (green triangles) , and sp i ral (bl ue stars) galaxies as a function of B-band luminosity, 
with selection weights applied. Binomial errors are computed using thelWilsonl (|1927k method. 

only the central SO frequency seems to be independent 
of halo mass. 

Halo mass also has a strong effect on satellite galaxies, 
although in a somewhat different fashion. As halo mass 
grows, the fraction of satellites which are spirals falls, just 
as happens for central galaxies. But satellite spirals in 
higher-mass halos are clearly being replaced by lenticular 
galaxies, which become the dominant type of satellites 
for M ha]o > HFh^Mm. 



fractions of central and satellite galaxies depend on halo 
mass, we see some striking differences. This is not actu- 
ally true for spiral galaxies, which decline in frequency as 
halo mass increases for both central and satellite galax- 
ies, just as we saw for galaxies in general (Figure [6]). 
But when we look at elliptical galaxies, we see a clear 
dichotomy: the elliptical fraction for satellite galaxies is 
roughly constant = -0.34 ± 0.93, P (3=0 = 0.72), while 
the fraction for central galaxies is a steeply increasing 
function of halo mass ((3 = 2.2±0.38, P p=0 ='8.2xI0~ 9 ). 

SO galaxies also show a dichotomy: the fraction of cen- 
tral galaxies which are SO is roughly constant (the appar- 
ent decrease for Mh a .i < 10 12 h _1 M m is li kely a result 



of the selection effect discussed in Section 2.8 above) 



but the fraction of satellite galaxies which are SO jumps 
from Oto 6 % to 69.0^% in higher-mass halos! This 
appears to be a highly significant difference; the logistic- 
regression value of P / s=o = 0.035 probably understates 
the significance because an abrupt transition like that 
seen here is not well modeled by the logistic curve (it 
would, of course, not be well modeled by a simple linear 
fit to the binned fractions, either). 

One clear implication from this analysis is the impor- 
tance of group halo mass for determining galaxy mor- 
phology. For central galaxies, halo mass is clearly more 
important than galaxy mass. The frequency of central 
spirals is high and essentially constant over the stellar 
mass range M* = 10 10 5 -10 12 M Q . But when we look at 
variations in halo mass, we can see central spirals being 
replaced by central ellipticals as the halo mass grows; 



lhalo > 

3.5. Co-evolution of Morphology, Star Formation and 
AGN Activity 

Physical processes which transform spiral galaxies into 
elliptical or lenticular galaxies should describe the de- 
struction of a galaxy's disk in the case of ellipticals, the 
enhancement of the bulge component for both ellipticals 
and bulge-dominated lenticulars, and the dissolution of 
spiral arms (lenticulars). Such processes should lead to 
the observed correlations between morphology and en- 
vironment. In addition, these processes may be respon- 
sible for circumnuclear starbursts and AGN, along with 
the removal of gas and the general suppression of star 
formation. 

We use the spectroscopic information from the SDSS 
survey to characterize the ongoing star formation and 
nuclear activity in our sample. Ha emission-line flux 
correlates strongly with the star formation rate, and is 
strongly bimodal for galaxies with and without signif- 
icant star formation or nuclear activity ( |Balogh et al. 
2004|. SDSS spectral fibers sample only the central 3" 
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11.8 



12.0 



-19 galaxies as a function of stellar mass, 



(diameter) of each galaxy, corresponding to 0.61-2.35 
kpc at z = 0.01-0.04. This undersamples the galaxy as a 
whole, and will miss much of the ongoing star formation 
and Ha emissi on at large radii. 



We use the |Brinchmann et al. ( 2004| hereafter B04) 
calibration of Ha equivalent width, which corrects the 
emission flux for underlying stellar absorption. We find 
a bimodality in this quantity, with a low emission peak at 
-0.5A ( emission is negative), which we presume should 
be at OA (no emission for a truly passive galaxy). We 
define a galaxy as "passive core" if the Ha equivalent 
width is within 2.5rr of this peak (where a in this case is 
the error in equivalent width, as estimated by B04, scaled 
up by a factor 2.473 calibrated to repeat measurements of 
the same galaxy r^j . The typical value of 2.5c is ~ 0.6A 
for these galaxies. 

The upper panel of Figure [8] shows the fraction of disk 
galaxies with a passive core as a function of Hubble type. 
Passive-core galaxies are, unsurprisingly, most common 
amongst SO galaxies, and the fraction decreases going 
from the RC3 type S0~ through SO to S0 + , going from 
E-likc to spiral-like SOs. Among spirals, Sa galaxies are 
much more likely to have passive cores than later Hubble 
types (~ 25% Sa, ~ 15% Sab, ~ 5% Sb, going to zero 
for later types). The anomalously low SO/a passive-core 
fraction is difficult to explain, and we avoid interpreting 
this for now. 

The lower panel of Figure [8] shows the fraction of SO/a- 

11 http:/ /www. mpa-garching.mpg.de/SDSS/DR4/raw .data. html 



Sb galaxies with a passive core as a function of the RC3 
value r25. This is the logarithmic axial ratio log 10 a/6, 
where a and b are the semi-major and semi-minor axes 
measured out to a surface brightness of 25 mag arcsec -2 
in the B-band. Low values of r25 imply near circular 
isophotes. The core of an early- type spiral is much more 
likely to be passive for a face-on inclination (low axis ra- 
tio), for which there will be no fiber contribution from 
the outer disk. In such early-type, face-on spirals the 
SDSS fiber spectrum can be dominated by bulge light. 
These correlations therefore suggest that star formation 
in the disk can be either truncated or heavily obscured 
in the inner regions, so that no significant Ha emission 
is detected in the fiber. Since we find a very low pas- 
sive fraction in highly inclined, early-type spirals (where 
the outer disk is likely to be projected into the SDSS 
fiber aperture), we infer that the outer disks of early- 
type spirals are typically still forming stars, and that 
most passive-core spirals are passive only in the inner re- 
gions. This is consistent with visual inspection of color 
JPEG images of these galaxies: in many cases the in- 
ner bulge+bar region appears red, with blue outer spiral 
features, often separated by a ring. 

The resulting picture of an inner truncation for star 
formation is a necessary simplification, given the limita- 
tions of our data. For comparison, longslit spectroscopy 
of four spiral galaxies lacking emission in the SDSS fiber 
revealed emission in the outer regions of two - and strong 
Balmer line ab sorption indicated re cent star formation in 
the other two pshigaki et al.||2007[ ). 
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Star formation is not the only possible source of Ha 
emission, of course. Ha emission traces ionized gas - 
and therefore requires both gas and ionizing radiation. 
The nature of that ionizing radiation can be explored via 
emission-line ratios, tracing the relative importance of 
different ionization levels and subsequent transitions. In 
Figure [9j we use the [N II]A6584/Ha vs [O III]A5007/H/3 
diagnostic line-ratio diagram, commonly used to sepa- 
rate normal star forming galaxies from harder radiation 
fields typical of Seyfert and LINER (low -ionization nu- 
clear emission-line region) type galaxies (Baldwin et al 



lines) > 3). Overplotted are galaxies in the RC3-SDSS 
sample with S/N(all lines) > 2, keyed b y morphology 
(see ca ption) . Anything to the right of the |Kewley et al. 
( |2001[ ) dashed line cann ot be explained by""horm al star 



1981[ the "BPT" diagram). The small points m the 
background trace the parent sample, demonstrating the 
overall distribution of galaxy line-ratios (where S/N(all 



burst models, whilst the Kauffmann et al. (2003) dashed 
line demarcates the boundary ot normal star forming 
galaxies (to the left). It should be noted that RC3-SDSS 
galaxies are at lower redshifts than the typical parent 
sample galaxy, and so fiber spectra will be more domi- 
nated by nuclear emission. Nonetheless, it is interesting 
that almost all E/S0 galaxies with emission are classi- 
fied as LINERs ([O III]/H/3 < 3 implies a softer- than- 
Seyfert ionization field). Some SOs extend to the region 
between the two dashed lines (composite/transition sys- 
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inconsistent with starbursts, but with [O IIlJ/ri/3 ratios typically lower than than ii, suggestive of LINERs rather than Seyfert galaxies. 
Recall that the parent sample goes to higher redshift, which means that SDSS fibers sample light out to larger physical radii. 



terns, pCaufFmann et al.|20"03 ), whilst SO/a-Sb (and some 
later-type) spirals extend from the bottom (high metal- 
licity) end of the star forming locus, up through the com- 
posite and LINER region to the Seyfert regime (with 
LINERs again the dominant population). As we sample 
only the central < 1 kpc, the ionization source may well 
be related to accretion onto a super-massive black hole 
(SMBH). 

We classify galaxies with emission lines into two basic 
categories. Anythi ng with broad emissio n lines or lying 
to the right of the Kewley et al. (20011 line is termed 
' AGN" . All other emission line galaxies are called "star- 
forming". Whilst these are commonly used definitions, 
we note that LINER-like ionization may also result from 
older stellar populations, sh ocks or interaction with hot, 



(Very Large Array) FIRST (Faint I mages of the Radi o 
Sky at Twenty-Centimeters) survey ( |Becker et al.|1994[ ). 
This 1.4GHz survey has a nominal detection threshold 
of 1 mJy, with a 90% confidence positional error circle of 
radius 1" (0.5" at the 3 mJy level). We first identified 
"nuclear" sources by requiring a match between SDSS 
and FIRST positions within 2"; this yielded 261 FIRST 
sources with fluxes > 1 mJy matched to the SDSS-RC3 
sample. 

Continuum radio emission at 1.4GHz originates via 
synch rotron emission, either in supernovae remnants 



Condon & Yin 1990 Condon 1992 Weiler et al 



2002 i - which correlates with star formation - or m m- 



X-ray emit ting gas (see e.g. Sarzi et al.||2010 Capetti & 
Baldi|20lT I and not solely from accretion onto a SMBH 



teractions between jets and the amb i ent medium (e.g. 
Burbidge|1956| |Guthmann et al.|2002") |Kaiser|2006[ ). We 



lb understand the possible role of AGN feedback for 
the suppression of star formation, it is also useful to ex- 
amine the radio properties of these galaxies. Models 
of galaxy formation invoke radio-mode AGN feedback, 
which suppresses cooling onto galaxies living at the cen- 
ters of massive halos by coupling t he kinetic energy of 



a radio jet to the cooling gas (e.g. Bower et al. 2006 



Croton et al. 2006). To examine the role of radio- mode 



feedback, we cross-correlated our sample with the VLA 



are primarily concerned with early-type galaxies, most 
of which are not forming stars. However, to ensure that 
we are dealing with nuclear radio sources which are un- 
likely to be due to star formation, we restricted ourselves 
to those galaxies whose radio emission was significantly 
stronger than what one would predict from the opti- 
cally determined star formation rate; we refer to these 
as "radio-AGN" sources. Specifically, we used the fiber- 
based B04 star-formation rates (which are insensitive to 
aperture corrections and well matched to the average 
compact nuclear radio source), and then estimated the 
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expected s tar-formation-based r adio luminosity using the 
relation of Hopkins et al. (2001): 



£i4GHzSF = 1-26 x 10 21 x SFR(M yr- 1 ) WHz~\ 

, . . . (3) 

converted to a Kroupa (2001) IMF for consistency with 
B04. Since Hopkins et al. found approximately an or- 
der of magnitude scatter in their relation, we impose a 
conservative limit of £i.4GHz, first > 10 x £i.4GHz,SF in 
order to identify bona- fide radio- AGN sources. As an 
example, a source at the detection limit of the FIRST 
survey (flux of 1 mJy) at the high-rcdshift end of our 
sample (z ~ 0.05) could be explained by a central SFR 
of ~ 7.5 M Q yr _1 , but would only be counted as a radio- 
AGN if the measured SFR was < 0.75 M yr" 1 . 

The resulting 140 galaxies in the SDSS-RC3 sample 
for which we found radio- AGN sources are indicated in 
Figure [9] by the black dots inside the galaxy symbols. 
Most elliptical and SO galaxies with radio- AGN sources in 
our sample are LINERs, which suggests an active SMBH 
is indeed pres ent. 

Figures 10 to 13 show the fractions of elliptical (red 



circles) and SO (green diamonds) galaxies which belong 
to each spectroscopic class (passive-core, star-forming, 
AGN), or with radio- AGN sources, divided into central 
(solid symbol) and satellite (open symbol) categories. 
This fraction is plotted both against stellar mass M* 
(left panels) and halo mass Mhalo (right panels). These 
figures illustrate how the nuclear spectroscopic and ra- 
dio properties of galaxies depend upon their morphology, 
mas s an d en viron ment. In contrast to the plots in Sec- 
tion 3.2 and |3.4[ we extend the stellar-mass plots here 
down to 10 M , becaus e th e bias towards spirals at 
low stellar masses (Section|2.8|and Figure[3l is no longer 



relevant. 

Figure [10] shows that the cores of elliptical galaxies 
(both central and satellite) are much more freqently pas- 
sive than SOs (70^g% for central Es compared to 34±5% 
for central SOs, 90tg% for satellite Es compared to 65tg% 
for satellite SOs). The fraction of passive-core SOs is 
highest for low mass satellites, whi ch a ll live in massive 
Mhalo > lO 13 h _1 M halos (section [O) . 

Figure [TT] shows that the only significant population 
of early-type galaxies with core spectra indicating star 
formation are central SOs with low stellar masses in low- 
mass halos (/? = — 9 ± 4.1, Pp=o = 0.029 versus M* and 
P = -8.6 ± 4.1, P^o = 0.04 versus M ha i ). Selection 
effects are relevant, such that low mass star-forming SOs 
only make it into our 5-selected sample because they 
are bright in that band. However, this population ex- 
tends into bins of stellar mass above our threshold M* > 
10 10 ' 5 Mq, suggesting a physical truncation of star forma- 
tion (and/or the presence of harder, AGN-like ionization) 
in central SO galaxies of halos Mhalo > lO 12 h _1 M . 
There are no notable trends for star-forming cores in el- 
liptical galaxies. 

Figure 12 shows that emission-line AGN (mostly LIN- 
ERs, Figure Km are found much more frequently in SOs 
than in ellipticals (combined central plus satellite pop- 
ulations: 9 ± 3% for ellipticals and 34 ± 4% for SOs), 
and their frequency increases with stellar mass (no- 
tably for SOs, although this is not formally significant: 
j3 = 2.8 ± 1.5, Pj3 = a — 0.067). However there is no mea- 



surable dependence on environment for AGN fraction, 
either in terms of halo mass or central versus satellite 
status. 

Figure [13] further explores the AGN theme in terms of 
radio emission. The fraction of galaxies with radio- AGN 
sources increases with both galaxy mass and halo mass. 
For central ellipticals, we find (3 = 2.8 ± 1.4^,3=0 = 
0.052 for radio-AGN fraction versus galaxy mass and 
p = 1.5 ± 0.53, P^ =0 = 0.0074 versus halo mass. The 
right-hand panel of the figure shows that for a given halo 
mass, central galaxies are more likely to host radio-AGN 
sources than satellite galaxies; the least likely hosts are 
satellite SOs. 

3.6. Dependence on Environment and Hubble Type of 
Bar and Ring Fractions 

Roughly two-thirds of local spiral galaxies - and a 



smaller fraction of SO galaxies - are barred (e 
et al.||2000[ |Menendez-Delmestre et al.)|2007[ ) 
A' -body simulations have long shown that ba 



Eskridge 
Although 

body simulations have long shown that bars can form 
spontaneously in isolated disks, simulations have also 
shown that bar formation can be triggered by tidal in 
g.,|Nogu chi||19871 |Salo||1991| |Noguchi||1996 
■ PUT" ' 



teractions (e. 



Berentzen et al. Tf is therefore plausible that loca 



environment might influence the frequency (and possibly 
the size or strength) of bars in disk galaxies. Similarly, 
although outer rings are well understood as being pri- 
marily due to the interaction of a bar' s Outer Lindblad 



Reson ance with gas in the disk (e.g., Buta & Combes 
|1996| ), the fact that they are features of the outer disk 
means they are in principle more vulnerable to inter- 
actions than other, more central structures. Thus, we 
might expect that local environment could also influence 
the frequency of outer rings. 

To estimate the bar fraction, we consider both strong 
(RC3 class SB) and weak (RC3 class SAB) bars in disk 
galaxies (SOs and spirals considered separately) . We also 
restrict the sample to relatively face-on galaxies: those 
with RC3 axis ratios a/b < 2.0 (r25 < 0.301). The latter 
restriction excludes highly inclined galaxies, where op- 
tical bar detection becomes difficult or impossible. We 
warn the reader in advance that we are probably underes- 
timating the true bar fraction, since some bars will have 
been misse d due to dust obscuration (see, e.g., Eskridge 
et al. |2000[ ) and weaker and smaller bars will have been 



difficult to identify due to resolution effects. Our analysis 
should thus be seen as investigating the possible effects 
of group environment on large, strong bars, rather than 
on all possible bars. 

We do find some evide nce for trends in bar fraction 
with galaxy mass. Figure 14 suggest that the spiral bar 
fraction increases with galaxy stellar mass, although the 
significance of this depends on which set of stellar masses 
we use (e.g., P/3=a = 0.020 using G05 masses, but only 
0.058 with our preferred Z09-based masses, and 0.047 
for the Y07 masses). A similar trend trend may exist 
for central SOs (Pp=v = 0.017 using Z09-based masses, 
with similar values for the other mass estimates ). Such 
a trend would be consistent with the findings of |Nair fc] 



Abraham (2010b) for a similar mass range. 

We find no evidence for trends in bar fraction versus 
halo properties, either as a whole or when considering 

right 



14 



central and satellite galaxies separately (Figure 
panel) . This is consistent with other recent studies of bar 
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Fig. 10. — Fractions of Mg < — 19 elliptical (red circles) and SO (green diamonds) galaxies with passive cores, separately for central 
galaxies (filled symbols) and satellite galaxies (open symbols), using selection and V/V m &x weights. This is presented in bins of stellar mass 
(left panel) and halo mass (right panel). Ellipticals, and low mass satellite SOs of high-mass halos, have predominantly passive cores. 
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fractio n with environment. For example, Marinova et al. 
( |2009[ ) found no difference in bar fraction (determined 
using optical HST images) across a range of lo cal envi 



ronments in the Abell 901/2 Supercluster, and Aguerri 
et al. (20091 found no evidence for a dependence ot bar 



traction (determined from SDSS images) with local envi 
ronment; both of these studies used local proje cted sur 
face d ensity of galaxies as the "environment" . |Li et al. 



( 2009 ) also found no difference in the clustering proper- 
ties of barred versus unbarred galaxies. 

For outer rings and pseudorings, we restrict ourselves 
to SO-Sbc galaxies, which is where almost all such rings 
are found (Buta & Combes 1996). We also consider 
rings and pse udorings separately, in p art because pre- 
vious work by Elmegreen et al. (1992) suggested there 
might be divergent trends tor the two subtypes. There 
is no clear evidence for any trend of outer ring or pseu- 
doring frequency with halo mass. The right-hand panel 
of Figure [15] appears to suggest a decreasing frequency 
with highernalo mass, at least for central galaxies, but 
this is not statistically significant, even if we lump outer 
rings and pseudorings together {Pp=o — 0.97 for central 
galaxies, Pp = o — 0.25 for all galaxies). What does seem 
to be present is a decrease in outer ring or pseudoring fre- 
quency with increasing galaxy stellar m ass, at least for 
central galaxies (left panel of Figure |15[ ). These trends 
are relatively shallow, but statistically significant (e.g., 
Pp=o = 0.00011 for central galaxies, lumping both outer 
rings and pseudorings together). Since outer rings ar e 
usually associated with bars (e.g., |Buta fc Combes|1 996), 
and since the frequency of bars apparently increases with 
stellar mass as noted above, the lack of outer rings in 
more massive galaxies is highly significant. 

The (tentative) absence of environ mental trends for 
outer rings appears to contradict what |Elmegreen et ah] 
( 1992 ) found: they argued that outer rings decreased in 



frequency for denser environments, while the frequency 
of pseudorings increased. One difference is that their 
main analysis was restricted to strongly barred SO + 
SO/a galaxies only, and it is not clear how to compare 
their "field", "pseudo-field" (possible group members), 
"group" , and "binary" classifications with our group halo 
properties. For example, a "binary galaxy" system could 
be low-halo-mass group with two significant members, or 
a subset of a larger group with a higher halo mass. It 
is also worth noting that the high outer-ring fraction re- 
ported by Elmegreen et al. for field galaxies is based on 
very small sample sizes (3 field galaxies, 4 pseudo-field 
galaxies) . 

The RC3 catalog also provides classifications for inner 
ring/spiral structure in spirals, specifying whether the 
main disk (outside the bar, if present) is purely spiral (s), 
contains an inner ring (r) , or has an intermediate broken- 
ring or pseudoring appearance (rs); see, e.g., Figure 1 of 
Buta et al. (1994). We have looked for possible correla- 
tions of r/rs/s frequency with environment in our sample, 
but find no evidence for any clear trends with hal o m ass, 
dispersion, or number of group members (Figures 16 and 
17). We note that the lack of any tren d for the mner- 
rmg fraction is possibly in conflict with Madore (1980), 
who found that galaxies with inner rings had tewer close 
companions (galaxies at projected distances < 50 kpc) 
than inner-spiral or inner-pseudoring galaxies. However, 
Madore pointed out that the number of close compan- 



ions, so defined, seemed independent of whether or not 
a given galaxy was in a group. 

We do find that the frequency of inner rings increases 
with galaxy luminosity and mass (/3 = 2.4±0.36, P^—q — 
1.8 x 10~ n for M*, with similar slopes and similarly 
small values of P/3 = o for G05 and Y07 stellar masses). 
Given t hat inner rings are more com mon in early-type 
spirals (de Vaucouleurs & Buta 1980), this trend could 
be a side effect of the tendency of early type spirals to 
be more massive than late-type spirals; it could also be 
due to the apparent increase in bar fraction with stel- 
lar mass noted above, sin ce inner rings are usual ly due 
to bar-related resonances (Buta & Combes 1996). Fig- 
ure [FT] does suggest a corresponding decrease m inner- 
spiral fraction with galaxy mass, though this is not sta- 
tistically significant (ft = -0.7 ± 0.61, P^ =0 = 0.25). 



4. DISCUSSION 



4.1. 



Where and How Are Elliptical Galaxies Formed? 

Figure [5] shows that the global fraction of elliptical 
galaxies increases with stellar mass, while Figure [6] shows 
no such trend with halo mass, velocity dispersion or num- 
ber of group members. This could naively be interpreted 
as evidence for a purely mass-dependent formation of 
ellipticals, independent of environment. However, Fig- 
ure [7] shows that a strong trend with halo mass does 
exist when only central galaxies (the most massive ones 
in their group) are considered. No such trend is seen 
for satellites; the fraction of ellipticals for satellites cor- 
relates only with stellar mass. The similar stellar mass 
dependences of central and satellite ellipticals, and the 
relatively infrequent occurance of satellite ellipticals, ar- 
gues for a common parent population for both types. 

Simulations show that disks are largely destroyed by 
major merge rs, leading to the formation of elliptical 
galaxies (e.g. Barnes 1988). These events preferentially 
occur at the centers ot halos: dynamical friction brings 
satellite galaxies to the bottom of the potential well 
where they merge with the central galaxy. Cosmologi- 
cal simulations show that mergers between the subhalo s 
hosting satellite galaxies are rare ( |Angulo et al. 2009). 
The highest mass halos have the richest halo merger his- 
tory, and the central galaxies of these halos are the most 
massive galaxies in the Universe (the most extreme ex- 
amples being cD galaxies in clusters), with the most ex- 



De Lucia & 



tensive merg er histories of any galaxy (e.g 
Blaizot|[2007 ). 

It all ellipticals are formed as (or transformed into) el- 
lipticals while they are still central galaxies within their 
own halos, then the observed correlations are perfectly 
consistent wi th our understanding of how ellipticals form 
via mergers ( |De Lucia et al.||2011[ ). The more massive a 
halo, the more likely its central galaxy is to have under- 
gone multiple major mergers, and thus the more likely 
it is to be an elliptical; such galaxies will naturally also 
tend to be more massive. Satellite ellipticals are galaxies 
which formed as ellipticals at the center of their pro- 
genitor halos - thus partaking in the general trends just 
outlined - and were subsequently accreted as ellipticals 
onto their current halos. Correlations with stellar mass 
will therefore persist for satellite ellipticals, whilst their 
previous (central-galaxy) correlation with halo mass is 
lost. 
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Fig. 14. — Fractions of strong + weak bars in (Mg < —19) SO (green triangles) and spiral galaxies (blue stars) as a function of the 
galaxies' stellar and halo masses, with selection and V/Vma,* weights applied. Filled symbols indicate central galaxies and hollow symbols 
indicate satellite galaxies. In this plot, we consider only galaxies with isophotal axis ratios a/b < 2 (i.e., relatively face-on). 
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Fig. 15. — As for Figure [14] but now showing fractions of SO-Sbc galaxies which contain either outer rings (black circles) or outer 
pseudorings (green squares) as a function of the galaxies' stellar and halo masses. Filled symbols indicate central galaxies and hollow 
symbols indicate satellite galaxies. 
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Fig. 16. — As for Figure [l4| but now showing fractions of spiral galaxies which contain inner rings (RC3 "r" classification, black circles) 
or inner pseudorings ("rs", green circles) as a function of the galaxies' stellar and halo masses. In this plot, we consider only galaxies with 
isophotal axis ratios a/b < 2 (i.e., relatively face-on). Filled symbols indicate central galaxies and hollow symbols indicate satellite galaxies. 
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In this scenario, most properties of elliptical galaxies 
(in addition to mass and morphology) are determined 
by their formation as central galaxies, and subsequently 
frozen at the time of their accretion onto larger halos as 
satellites. Thus, aside from possible differences in mean 
stellar age, satellite and central ellipticals should follow 
the same general scaling relations. Studies have found 
that the well-known elliptical-galaxy scaling relations are 
indeed largely independent of environment, These in- 
clude the slope of the fundamental plane (e.g.,|de la Rosa| 
let al. 1120 011 IReda et al.|[2005l iBernardi et al.||2TO| bu t 
TrOnofrio et al.|l2008( |La Barbera et al.||2010D, 



see also 



the color-m agnitude and Rorm endy relations (e.g!, Hogg 



et al. 2004; Rcda et al. 20 05]) , and luminosi ty-size rela- 



tions Wjm k al.|2010p . |guo et al.| ( |2009[ ) and 



Weinmann 

et al. 1 2009 ) studied galaxies using the Y07 group cat- 
alog, specifically contrasting mass-matched central and 
satellite galaxies, and found no size or structural dif- 
ferences for "early- type" galaxies. Note that many of 
these studies lumped elliptical and SO galaxies together, 
so there is in principle the possibility of confusion if el- 
lipticals and SOs have trends that happen to cancel out 
when they are combined. 

While there is evidence for age differen ces between el- 
lipticals in different environments (e.g., |Thomas et ah 
20051 |Bernardi et~aI1|2006| |Saglia et al | |2010p , thisls 
still consistent with the overall picture. Ellipticals in 
high density environments are more likely to be satellite 
galaxies within massive halos, with properties frozen at 
the time of the accretion. 

The picture outlined above implicitly assumes that 
once an elliptical has formed at the center of its halo, 
it is able to remain an elliptical. In the case of ellip- 
ticals created by "wet" (gas-rich) mergers, the problem 
is how to prevent significant residual gas from continu- 
ing to form st ars; quasar-mode feedback is a promising 



solution (e.g., Granato et al. 2004 Springel et al. 2005 
Hopkins et al. 2006| ). A more general, long-term prob- 
lem is that posed by the presence of hot gas in the halo. 
Since the center of a halo is the natural destination for 
halo gas that is able to cool, some mechanism must exist 
for suppressing such cooling - otherwise, cool gas would 
accumulate in the center of the halo and potentially form 
a new stellar disk. (Note that satellite ellipticals do not 
suffer from this problem, since they do not sit at the 
centers of halos.) 

This cool, low entropy gas has been observed in some 
massive clusters - but it is not ubiquit ous; clusters with 



no cooling flows ar e in the majority (e.g. Nesci 1991 Cav- 
agnolo et al.|2009 ) . The existence of cooling gas appears 



to be a requirement for ce ntral cluster galaxy to host 
either star formation (e.g. Donahue et al.||2010| |Hicks 



et al.||2010[) or radio AGN (e.g. |Sun||20^y TOras been 



proposed that energy from a radio jet can offset cool- 
ing in a cluster and that this suppresses the growth of 
the central galaxy. Cavities containing low-density, hot 
gas have been observed spatially coincident to the ra- 
dio lobes, and the total energy required to create such 
cavities is typically enough to balance cooling (e.g 
Namara et aTJl20"00l |Fabian et al.||20"00"l |Dunn fc Fa' 
2006| |Cavagnojo et al ||2010[ ) . 

If radio-mode AGIN feedback is also applicable for 
galaxies in lower-mass halos, then it could prevent sig- 
nificant cooling onto ~ galaxies, suppressing their 



Mc- 



bian 



growth and allowing models of galaxy formation to match 
the high- mass exponential cutoff seen in the ga laxy mass 
function ( |Bower et al.|2006 Croton et al.|2006 ). Cavities 
have been observed in the hot gas component of galaxy 
groups ( |Dong et al.||2010l and even individual elliptical 
galaxies (fBaldi ei al. 2009}. The high fraction of radio 
sources we find for central ellipticals and SOs in halos 



13 



see also 



with masses of 10 13 -10 14 - 5 h' 1 My ( Figure 
Best et al.|[2007||Pasquali et al.|2009"| ) is potentially fur 
ther support for the possibility that radio-mode feedback 
operates in halos of these masses. It is also possible that 
other forms of AGN feedback can suppress star formation 
in lower-mass halos; for example, Schawinski et al. ( 2009 \ 
have found evidence that low- luminosity AGIN activity is 
associated with the disappearance of central molecular 
gas in SO and elliptical galaxies. 

More detailed models will be necessary to determine 
whether the fractions of central galaxies with elliptical 
as opposed to disk morphology, and their dependence on 
halo mass, can be quantitatively explained in the context 
of the expected merger history and suppression of disk 
re- formation (Wilman et al., in prep). 

4.2. Where Are SOs Formed? 

As with ellipticals, the dependence on halo mass of 
the SO fraction is very different for satellite and cen- 
tral galaxies. The remarkable change in satellite SO fre- 
quency, which jumps from 01q' 6 % in lower-mass halos to 

69. 01^% for halos with masses > 10 13 /i _1 M , suggests 
that spiral arms are often suppressed in a disk galaxy 
once it is accreted onto a halo, if the halo is more mas- 
sive than 10 13 h- 1 M . 

In fact, we can argue that the majority of all present- 
day satellite SO galaxies became SOs after they were 
accreted into halos more massive than ~ 10 13 h~ x Mq. 
Since z = halos with M halo < IO^Ii^Mq do not 
have satellite SOs, any pre-existing SOs which formed in 
less massive progenitor halos and then fell into massive 
halos (Mhalo > 10 13 h _1 M Q ) must have originally been 
central galaxies. Assuming that these progenitor halos 
had a distribution of central morphological types similar 
to what we see today (if anything, they probably had 
lower SO fractions in the past), then the mean central 
SO fraction for the accreting progenitor halos would be 
~ /|™ = 20.1 ± 1.6%. The fraction of SOs in massive 
halos which are post-processed is therefore: 



/•post rsat 

Jso ~ Jso 



ft 



con 
SO > 



(4) 



where /gQ st is the fraction of massive-halo satellite galax- 
ies which are post-processed SOs, and /|q* is the fraction 
of massive-halo satellites which are SOs. The present-day 
SO satellite fraction (/go* — 69 ± 4%) thus requires that 
almost three quarters of these galaxies (i.e., 49 ± 5% of 
Mhalo > 10 13 h^ 1 M satellite galaxies, or 71 ± 8% of 
satellite SOs) fell into their present-day halos as spirals, 
becoming SOs during or after the accretion process. 

We can go one step further, and apply equation [3] as a 
function of galaxy stellar mass. We limit this exercise to 
the stellar mass range 10 10 5 M Q < M* < 10 115 M Q for 
two reasons: first, as previously noted, SOs are lost due 
to selection effects below 10 10 5 M Q ; second, there are no 
satellite SOs in our sample with masses > 10 115 M Q (see 
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Figure [7]). For simplicity, we assume that there is no sig- 
nificant stellar-mass change during or after the accretion 
process. The upper panel of figure [18|shows how /g™ and 
fso depend on stellar mass. We also show the equivalent 



and fig. 



What 



fractions for spiral galaxies: i.e., /g™ 
is striking about this plot is the apparent stellar-mass 
trend. For M* < 10 n M Q , the satellite SO fraction is 
clearly too high to be explained by the pre-processed SO 
population alone, requiring substantial conversion of ac- 
creted spirals into post-processed SOs. At the high-mass 
end, on the other hand, most or all of the satellite SOs 
can be explained as pre-processed SOs, with little or no 
conversion of spirals required. This is consistent with the 
fact that the satellite spiral fraction in the highest-mass 
bin is basically the same as the central spiral fraction, 
suggesting that ~ all of the highest-mass spirals have 
rema ined spirals after accre t ion. T his can be compared 
with van den Bosch et al. (2008 esp. their Figure 8), 
who apply a similar argument to suggest that the high 
fraction of massive, M* > 10 11 h -2 M galaxies which 
have red colors is largely due to their pre-processing as 
central galaxies. 

The lower panel of figure [18] makes some of this more 
explicit. Here we show our estimate for the fraction of 
pre-processed SOs, assumed to be equal to the fraction 
of central SOs in the top panel. We also estimate the 
fraction of post-processed SOs using equation [4j as hinted 
in the top panel, the fraction of satellite SOs which are 
post-processed increases to lower masses. Finally, we also 
plot the fraction of "missing" satellite spirals, which is 
the fraction of central spirals minus the fraction of satel- 
lite spirals in Mhaio > lO 13 h _1 M halos. This quan- 
tity represents the fraction of accreted spirals which are 
no longer present as spirals; these are assumed to have 
transformed into SOs or else merged with other galax- 
ies. In each mass bin, the fractions of post-processed SOs 
and missing spirals are equal within the errors, strongly 
suggesting that the missing spirals have indeed been con- 
verted into (post-processed) SOs. 



fore requires that star formation is suppressed, both 
by the exhaustion or removal of existing disk gas and 
by the prevention of further gas accretion from the en- 
vironment. Our results suggest a difference between 
the transformation of post-processed satellite galaxies in 
Mhaio > 10 13 h _1 M Q halos on the one hand, and that of 
central galaxies in halos down to our limiting halo mass 
of ~ 10 12 hr 1 Mq - independent of halo mass - on the 
other hand. We therefore expect different mechanisms 
to be responsible for SO formation in these two regimes. 

Finally, whatever mechanisms are operating, they also 
need to explain the observed structural differences be- 
tween spirals and SOs. SOs have traditionally been char- 
acterized as h aving extremely high B/T ratio s (e.g., 



Simien 
57 tor 



Dressier 198 0); for exam ple, the compilation of 
fc de Vaucouleurs] ( |1986[ ) has a mean B/T ~ ( 
SOs. More recent studies which account for variable 
bulge profiles (i.e., Sersic r 1 /" instead of de Vaucouleurs 
r 1 / 4 profiles) and the effects of additional components 
such as bars have resulted in lower B/T values for SOs 
- but the l arge study using 2D dec ompositions of disk 
galaxies by |Laurikainen et ah] ( |2010[ ) still finds that the 
mean (and maximum) B/T values for SOs are higher than 
those of spirals (though there is considerable overlap, and 
some SOs have B/T < 1); se e Figure 4 of that paper. 



Christlein & Zabludoff (20041 constructed separate lu- 
mmosity functions tor disks and for bulges as a function 
of B/T. For intermediate, 0.2 < B/T < 0.7 galaxies 
(assumed to be typical for SOs), they found that while 
the characteristic luminosity of bulges L* bulge strongly in- 
creases with B /T, the characteristic luminosity of disks 
L* disk is almost constant. Simple toy models are then 
used to show that this is inconsistent with changes in 
B/T resulting from pure disk fading - but that a model 
in which B/T evolves through the gro wth of bulges is 
consis tent with data. We also note that IBurstein et al.l 
( |2005| } presented evidence indicating the if-band lumi- 
nosities of local SOs were, on the whole, too high for all 
of them to be explained by fading of gas-stripped spirals. 



4.3. How Are SOs Formed? 

The persistence and visibility of spiral arms clearly cor- 
relates with the presence of gas and star formation: spiral 
galaxy disks are forming stars, while SO disks are typi- 
cally passive (e.g., Figure [8]). The absence of spiral arms 
in SO galaxies can be explained by the combination of 
increased random motions of disk stars with age, which 
erases existing spiral patterns, and lack of young stars in 
regular, "cold" orbits, which would otherwise maintain 
or ref orm spiral patterns (see e.g. section 6.1 of Sellwood 
2011| and references therein). 

Star forming galaxies at z ~ typically have enough 
atomic and molecular gas to maintain star formation 
for only another ~ 3 Gyr on average at present rates 
(for t he statistics of local galaxies see Saintonge et al. 
2011| . Ongoing disk star formation for times of order 
the Hubble time thus requires the availability of addi- 
tional gas, accreted from the surroundings. This gas is 
usually assumed to be shock-heated upon accretion onto 
> 10 12 ft" 1 M Q halos, resulting in a "hot atmosphere" 
which subsequently (i n the absence of furt her heating) 
cools onto the galaxy ( White fc Frenk|[l991 ). 

The transformation ot spiral galaxies into SOs there- 



4.3.1. Satellite SOs 

We have presented evidence that the majority of SOs 
in Mhaio > 10 13 h" 1 M Q halos were likely accreted as spi- 
ral galaxies, and have since been post-processed, leading 
to their current SO morphology. We have also shown 
that lower mass SOs (M* < 10 11 M Q ) are more likely 
to have experienced post-processing, whilst the fraction 
of higher mass SOs is more consistent with an accreted 
field population (pre-processed into SOs as central galax- 
ies). Thus, post-processing is more relevant for lower 
mass disk galaxies, while higher-mass accreted spirals are 
more likely to persist as spirals in the group environment. 
This puts potential limits on the post-processing mecha- 
nism^); in this section, we focus on those which remove 
gas from the galaxy, since that is an essential prerequisite 
for preventing further star formation and the persistence 
of spiral structure in the disk. 

A promising way of removing gas from spiral galaxies 
in massive groups and clusters is to remove it via interac- 
tions with the gas of the hot intra-group or intra-cluster 
medium (IGM/ICM). The physical nature of this interac- 
tion can - very broadly - take three forms: ram- pressure 
stripping ( Gunn fc Gott|1972 ) , viscous stripping ( |Nulsen| 



RC3-SDSS Galaxy Morphology and Environment 



23 




A pre-processed SOs 
O post-processed SOs 

missing satellite spirals 




10.8 11.0 11.2 

Galaxy Stellar Mass (log 10 M ) 



Fig. 18. — Top: Observed fractions of SO and spiral galaxies as a function of stellar mass M* for central galaxies (solid points and 
lines) and for the satellites of massive ( Mhaio > 10 13 h — 1 M©) halos (open points and dashed lines). Bottom: Estimated fractions for 
the satellite population of massive halos, divided into "pre-processed" SO galaxies (pre-existing SO galaxies accreted by the halo) and 
"post-processed" SO galaxies (accreted as spiral galaxies and subsequently transformed). These fractions assume that the morphological 
fractions of galaxies at the time of accretion matches the currently observed cen4rai-galaxy fractions (top panel). We also show the fraction 
of "missing satellite spirals" — these are galaxies accreted as spirals which are no longer present as spirals (i.e., they are presumably either 
transformed into SO's or merged to form ellipticals). 
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1982 1 , and thermal evaporation ( |Cowie &; Songaila" 1977 ) . 
Ram-pressure refers to the pressure exerted due to the 
galaxy's motion relative to the IGM/ICM and acts very 
quickly in regions of very den se, hot gas such as th e cores 



of massive clusters (see e.g. Moran et al. 2007); ram 
pressure stripping can potentially be strong enough to 
remove tightly-bound cold gas in the galaxy disk. Vis- 
cous stripping refers to the slower removal of low density 
gas via turbulence at the interface between the galaxy 
and IGM/ICM. Both of these mechanisms depend pri- 
marily on the IGM/ICM density and the galaxy veloc- 
ity. Evaporation instead removes gas via thermally in- 
duced collisions, and depends on the temperature of the 
IGM/ICM; it can also act on central galaxies, and so is 
not exclusive to satellite galaxies. 

There is direct observational evidence for ram-pressure 
stripping of col d gas from disk galaxies in the Virgo 
Cluster (see e.g. |Gavazzi et al.| 2008; Ch ung et al.||2009| 
and references therein) and also m the Coma Cluster 



(e.g., Vollmer et al. 20011. However, Virgo and Coma 
are relatively massive systems (clusters with Mhalo ^ 
10 14 1i _1 Mq), supporting dense ICMs and rapid galaxy 
motions, both of which lead to stronger galaxy-ICM in- 
teractions. Even here, stripping of disk gas in most cases 
seems to be partial, acting only on the outer parts of 
the disk and leaving a more tightly bound core gas com- 
ponent, which can cont inue forming stars. This is also 
seen in simulations (e.g. |Kapferer et al.|2009 1 . The main 
problem is that it is not clear whether stripping of tightly 
bound cold gas can ever operate effectively in lower-mass 
clusters and g roups (e.g., M halo ~ 10 13 -10 14 hr 1 M ). 
Moran et al. (2007) studied passive spirals and SOs in 
two clusters at z ~ 0.5 and concluded that ram-pressure 
stripping was significant only in the more massive clus- 
ter, with its much denser ICM. 

A more widely ap plicable process is "strangulation" 
(Larson et al. 1980[ ), which refers to the ram-pressure 
removal of just the hot gas halo of a galaxy. Since the 
halo gas is much thinner, hotter, and much less strongly 
bound than the cold disk gas, it is vulnerable to removal 
by lower velocities and lower ICM/IGM densities, mak- 
ing it a plausible mechanism for lower-mass clusters and 
groups. Simulations of strangulation find that, as for 
ram-pressure stripping of the cold gas, the central, most 
strongly bou nd hot gas component can survive the strip- 
ping process jMcCarthy et al.|2008||Kawata & Mulchaey| 
2008 |Bekki|[2(J09[ ). This is not necessarily a problem, 
however. Partial stripping of the hot gas causes star for- 
mation to be suppressed over longer timescales than if 
all of the hot gas and most or all of the cold gas were 
stripped (e.g., in the standard ram-pressure-stripping 
scenario). With little or no cooling of hot gas onto a 
satellite galaxy, the remaining cold disk gas will be ex- 
hausted in ~ 3 Gyr. 

Alternatively, a gravitational mechanism might be 
responsible for the post-processing of group galaxies. 
Galaxies with an extended history of minor mergers are 
likel y to form a high B/T remnant with SO morphology 



(e.g. Bournaud et al. 



2007). However, if galaxies follow 



dark matt er, the~p robability of satellite-satellite merg- 



ers is low (jAngulo et aLp009|) and so we expect this to 
th 



be more important in the case of central SOs (see be- 
low). A more promising, slow-acting mechanism is the 



cumulative effect of low-v elocity tidal interaction s with 
the other group galaxies (Bekki & Couch 2011). This 
is distinct from the frequent, high-velocity encounters in 
more massive clus ters, o ften referred to as "harassment" 
(e. g. |Moore et al ||1998l |Gnedin | |2003| ) . The simulations 
of |Bekki & (Jouch| (|201l| ) show that repeated interac- 
tions can drive tidal stripping and compression of the 
gas, leading to its more rapid exhaustion than in an iso- 
lated case; this process also produces bulge growth via 
central star formation - albeit by only A(B/T) < 10%. 
Tidal interactions also enhance the random motions of 
stars, which thickens the disk and helps suppress spi- 
ral arms. This mechanism acts preferentially on lower 
mass galaxies (which are more easily perturbed) and on 
galaxies in lower mass groups (M halo ~ 2 x 10 13 M Q ). 
Encouragingly, this is consistent with our population of 
post-processed SOs (predominantly M* < 10 11 M lower 
panel of Figure 18). 



4.3.2. Central SOs 



In most cases, the central galaxy of a halo lives at 
the bottom of the global potential well, and has little or 
no velocity offset with respect to the hot gas (see Skibba 
et al. 2011 for more detail and caveats). Therefore, strip 



ping cannot operate under these conditions, and we must 
look for other mechanisms to suppress star formation in 
central SOs. 

The centers of halos are where galaxy mergers predom- 
inate. Since mergers where one or more galaxy is gas-rich 
can induce starbursts that rapidly consume the gas (plus 
potential quasar- mode feedback) , this is a potential route 
for forming SOs. Minor mergers, in particular, should be 
more common than major mergers and have the addi- 
tional utility of tending to add mass to the bulge without 
completely destroying the di sk (Bckki 1998 Eliche-Moral 
et al.|2006||Bournaud et al.||2005[pl)O7| ; this would help 
increase the B/T ratio of central SOs. Of cours e, ma- 
jor mergers — a nd multiple minor mergers (e.g., Bour- 
naud et aL||2007 1 — are more likely to produce an ellip- 



tical remnant. Since higher-mass halos, with their richer 
merger history, are probably more likely to have had 
major mergers fo r their central galaxies (e.g., Wang & 



Kauffmann 2008), we would expect the central galaxies 
ot higher- mass halos to be ellipticals more often, which 
is indeed the trend we see for our sample (Figure [7]) . 

This suggests that the existence of SO galaxies with 
if-band lumin osities brighter than a ny spiral galaxy, as 



pointed out by Burstein et al. (2005), does not mean SOs 
cannot form through disk fading ol spirals. Instead, it 
suggests that the most massive SOs are a mixture of SOs 
which are still the central galaxies within their groups 
and formerly central (pre-processed) SOs which were ac- 
creted into massive groups and clusters. Less massive 
SOs are then more likely to be post-processed spirals. 

We are still left with the problem — similar to that we 
faced with explaining the lack of recent star-formation 
activity in central ellipticals — of how to prevent cool- 
ing halo gas from accreting onto central SOs and trig- 
gering new star formation. In principle, the same AGN- 
feedbac k me chanisms invoked for central ellipticals (see 
Section 4.1) could apply. We do in fact find relatively 
high radio- AGN frequencies in central SO's, especially 
in more massive halos (Figure 13). A possible problem 
is the lower masses of central black holes in SOs, which 



RC3-SDSS Galaxy Morphology and Environment 



25 



could make feedback less efficient. (Since SMBH mass 
scales with bulge mass rather than total galaxy mass — 
e.g., Kormendy & Gebhardt 2001, Kormendy et al. 2011 
— SOs will tend to have smaller SMBHs than ellipticals 
of the same ste l lar ma ss.) However, recent simulations by 
Gaspari et al. (2011) suggest that relatively weak feed- 
back may be all that is needed to keep cooling flows 
from developing in groups. We also see high fractions 
of optical AGN in central SOs as well (Figure 12 ), which 
could be helpful given the evidence that low-luminosity 
AGN may be associated with the disappearance of cen- 
tral molecular gas e ven in galaxies haying s tellar masses 
of a few xlO lo M ( |Schawinski et al.||2009[ ). 



the increase in SO fraction observed for clusters over the 



4.4. Evolution of the SO Fraction 

Comparisons with samples at different redshifts require 
care - especially where different selection limits, environ- 
mental definitions, and classifi cati on methods are em- 
ployed. The top panel of Figure [18] indicates that the SO 
fraction is sensitive to the mass or luminosity threshold 
imposed - a problem that is accentuated when compar- 
ing samples from different redshifts and with different 
photometry. 

Since most high-redshift studies where the SO fractions 
are computed are for massive clusters, and thus not good 
direct comparisons with our sample, we compare our lo- 
cal SO fracti ons with the z ~ 0.4 group and field sample 
presented by Wilman et al. ( 2009 ) . In Figure 5 of that pa- 
per, morphological fractions are computed with a lumi- 
nosity cut of My = —20.53, for comparison with cluster 
fractions from the literature^] Assuming passive evolu- 
tion, the luminosities of ellipticals and SOs will decrease 
between z ~ 0.4 and now; thus, we need to adjust our lo- 
cal luminosity cutoff accordingly. Based on t he updated 
calculations of van Dokkum fc Franx| (feOOlp j we esti- 
mate ~ 0.43 mag fading in V for E arid SO galaxies, and 
therefore adopt a local luminosity limit of My — —20.1 
for comparison purposesrj Note that star-forming spiral 
galaxies are not expectedto fade as much as elliptical or 
SO galaxies; some may even increase in luminosity. Con- 
sequently, the local sample may have an excess of spirals 
relative to the higher-redshift sample, since some spirals 
with My > -20.53 at z ~ 0.4 now have My < -20.1. 

The |Wilman et ~al~] p009| statistics are based on 
two overall classifications of environment: "groups" and 
"field". Their groups do not have halo mass esti- 
mates, but they do have velocity dispersions. Since 
almost all the Wilman et al. groups have dispersion 
> 175kms _1 , we divide our local sample into equivalent 
subsets, with "groups" defined as halos having disper- 
sions > 175 km s -1 (this amounts to 90 groups containing 
a total of 185 classified galaxies) and the "field" defined 
as all other halos. 

The SO fraction for the field shows no evolution over 
this redshift range: 10±|% at z ~ 0.4 versus 11.1+^% 
at z ~ 0.02. But in groups we do see some evidence for 
evolution: the fraction rises from 28 ± 4% at z ~ 0.4 
to 52.3^4Q% locally. This parallels at least qualitatively 

12 Selection and systematics are discussed in Section 4.2 of that 
paper. 

13 Available at | http://www.astro.yale.edu/dokkum/evocalc/| 

14 We compute V = g - U.mti(u - g) + 0.0885 plus the 
correction for oversubtracted background - see Section|2.5| 



same redshift range (e.g., |Dressler et al. 


1997 


Fasano 


et al.|2000 


Poggianti et al. 2009). Quantitatively, it also 



see ms to ag ree with the evolution in SO fraction observed 
by Just et al. (2010) for the lower-mass clusters in their 
sample (those with velocity dispersion < 750 km s" 1 ). 

4.5. Differences in AGN '/LINER Fraction Between 
Ellipticals and SOs 

Figure [12] shows that the fraction of emission-line 
AGN (mostly LINERS) is similar for central and satellite 
galaxies, and does not depend significantly on halo mass. 
It does increase with stellar mass, and is much higher in 
SOs than in elliptical galaxies of the same mass (except 
at the very highest-mass end). 

The mass of central supermassive black holes is known 
to increase with galaxy (more properly, bulge) mass (e.g., 
[ Magorrian et al.||T998| |Marconi fc BuiitpMl |Haring fe| 
|Rix|2004[ ). The increasing fraction of emission- line AGN 
with galaxy mass could thus potentially result from more 
powerful ionizing sources due to higher black hole masses. 
However, the mass of an elliptical galaxy is larger than 
the bulge mass of an SO of equivalent stellar mass. The 
implication of a tight bulge mass-black hole mass rela- 
tion is that an elliptical galaxy should also have a more 
massive black hole. In fact, however, we generally see a 
higher fraction of both optical AGN and radio- AGN in 
SO galaxies, (Figures 



12 and 



13) 



One possible solution might be that some SOs have 
retained a small amount of tightly bound gas in their 
cores, which can then be ionized by accretion onto the 
black hole. The absence in radio-bright elliptical galax- 
ies of Ha-emitting ionized gas suggests this gas is truly 
absent - either because it is too hot or because it has 
been expelled. Gas is heated to X-ray emitting te mper 



atures duri ng simulations of major mergers (e.g., Cox 
et al.||2006|) - but the inclusion of hot gas in s uites of 



merger sim ulations is still preliminary (see, e.g., Moster 
et al.||2011|), and the results are highly sensitive to the 



balance of heating and cooling. 

5. CONCLUSIONS 

We have created a catalog of 1064 nearby (median 
z ~ 0.02), bright (B < 16 and M B < -19) galaxies, 
which combines RC3 morphological classifications and 
the NYU-VAGC version of the SDSS DR4 photometric 
and spectroscopic data; magnitudes have been corrected 
to account for the undersubtracted background, and stel- 
lar masses were determined using q — i color-based mass- 



2009| ) and ap- 
he main mor- 



to-light ratios calibrated by Zibetti et al. 
plied to the i-band absolute magnitudes, 
phological classifications (elliptical vs. SO vs. spiral) were 
checked by visual examination of SDSS images; a total 
of 165 galaxies ended up being re-classified, and 55 more 
were classified for the first time. 

To this dataset we added halo masses and central vs. 
satellite status (where "central" = m ost massive galaxy 
in its group) from the group catalog of | Yang et al. ( 2007 



Y07), which resulted in a total of 911 galaxies with halo 
mass assignments; 729 of these are central galaxies and 
182 are satellites. The main advantage of the Y07 cat- 
alog is that it spans the full range of halo masses down 
to M ha i = 10 11 - 63 h~ 1 M Q . This allows us to describe 
the dependence of morphological fractions on halo mass 
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at z ~ for the first time, from smaller clusters down to 
single-galaxy halos. (We also determined total number 
of galaxies per group and group velocity dispersions for 
the Y07 groups in our sample, in order to check the ro- 
bustness of the halo-mass-based results.) Using the full 
SDSS DR4 catalog as a parent sample, we characterized 
the selection function of our catalog as a function of B 
magnitude (synthesized using SDSS colors) and galaxy 
size (rgo, the radius containing 90% of the Petrosian flux 
in r-band). Our sample is robust to a luminosity limit 
of Mg = —19.0, which corresponds to a stellar mass of 
~ 10 10 - 5 M for E and SO galaxies. Galaxies are 
weighted to correct for the selection bias, and are also 
weighted by V/V max to correct for Malmquist bias. This 
allows us to examine the fraction of galaxies of various 
types as a function of stellar and halo mass. We use 
a weighted logistic regression method to allow us to as- 
sess the statistical significance of apparent trends. This 
method has the advantages of modeling a binomial prop- 
erty without binning of data, and it allows for individual 
data points to have weights. 

We find that the global fraction of elliptical galaxies 
increases with galaxy luminosity and with stellar mass, 
but not with halo mass. The fraction of SOs declines to 
high stellar mass, but increases with halo mass, group 
velocity dispersion and number of neighbours. These re- 
sults are consistent with previous work at higher redshift 
which found the SO fraction to increase in groups rela- 
tive to the field populatio n, whilst the elliptica l fraction 
remains roughly constant (|Wilman et alj[2009 



The fraction of central galaxies with elliptical morphol- 
ogy increases with stellar and halo mass, consistent with 
their formation in mergers. In contrast, the fraction of 
satellite ellipticals is globally low at all halo masses, but 
increases with stellar mass. We interpret this as evidence 
that satellite ellipticals were formed as the central galax- 
ies of progenitor halos, which were subsequently accreted 
onto their present halo. 

Limited to M* > 10 10 5 M©, a modest fraction of cen- 
tral galaxies are SOs (20.1 ± 1.6%, with little dependence 
on stellar or halo mass). The remaining SOs are satellites 
of massive halos only - we find that the fraction of satel- 
lites with SO morphology rises from 0^q' 6 % in halos with 

Mhaio < lO^h- 1 M to 69.0±4 6% above this threshold. 
Presuming SOs to be spirals in which star formation has 
been suppressed (leading in turn to the suppression of 
spiral arms), we interpret our result as a strong indica- 
tion that there are two populations of SOs, in which star 
formation has been suppressed in different ways. 

Central SOs may be formed via suppression of star for- 
mation during minor mergers and/or by feedback from 
AGN, with similarities to the elliptical population. Satel- 
lite SOs which became SOs while they were still cen- 
tral galaxies within their progenitor halos constitute a 
pre-processed population which can account for up to 
~ 20% of the satellite SOs, including all of those with 
M* > 10 n M Q . 

However, the higher fraction of SOs in higher-mass ha- 
los implies that many satellite SOs were accreted as spiral 
galaxies. These accreted spirals were then post-processed, 
becoming satellite SOs, and are the dominant source of 
satellite SOs in the range 10 10 - 5 M© < < lO n M . 
Altogether, we estimate that 64 ± 11% of our satellite 
SOs were accreted as spirals. 



Central SO and elliptical galaxies frequently host radio 
sources, consistent with radio-mode heating of the sur- 
rounding hot gas. This heating may offset cooling onto 
these galaxies, and thus suppress star formation. How- 
ever, SOs host ionized gas components - mostly LINERs 
- much more frequently than elliptical galaxies of the 
same mass. 

We find no strong dependence of structural subcompo- 
nents - bars, inner rings/spirals, outer rings - on environ- 



ment, in con trast to some earlier studies (e.g. Elmegreen 



et al.||1992[ ), though we do find evidence that the ire 
quency ot both bars and inner rings increases, and the 
frequency of outer rings decreases, with galaxy mass 



By compa ring our results with the study of |Wilman 
et al. ( 2009 ) , we find tentative evidence that the fraction 
of bright SOs in intermediate-mass groups (those with 
velocity dispersions 175 km s -1 < a < 500 km s -1 ) has 
increased in the last ~ 4 Gyr, rising from 28 ±4% at z ~ 
0.4 to 52.3^4Q%. This is at least qualitatively consistent 
with increases in the fraction of SOs in clusters reported 
by other studies. 
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